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Abstract 

Thispaper reports the result of an experimentalstudyinthe 
NASA Glenn Icing Research Tunnel (IRT) to evaluate 
how well the MVD - independent effect identified previ- 
ously might apply to SLD conditions in rime icing 
situations. Models were NACA 0012 wing sections 
with chords of 53.3 and 91.4 cm. Tests were conducted 
with a nominal airspeed of 77 m/s (150 kt) and a num- 
ber ofMPD’s ranging from 15 to 100 pm with LWC of 
0.5 to 1 g/m 3 . In the present study, ice shapes recorded 
from past studies and recent results at SLD and Appen- 
dix-C conditions are reviewed to show that droplet di- 
ameter is not important to rime ice shape for MVD of 
30 pm or larger, but for less than 30 pm drop sizes a 
rime ice shape transition from convex to wedge to 
spearhead type ice shape is observed. Discussion on 
what may cause such transition is given and some evi- 
dence presented in this study suggests that the shape 
transition could be governed either by the accumulation 
parameter A c alone or by p 0 and A c together. 

Nomenclature 

A c Accumulation parameter, dimensionless 

AOA Angle of attack, ° 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 

c p Specific heat of air, cal/g K 

c P ,ws Specific heat of water at the surface tempera- 
ture, cal/g K 

d Cylinder diameter or twice the leading-edge 

radius of airfoil, cm 

h c Convective heat-transfer coefficient, 

cal/s m 2 K 

ho Gas-phase mass-transfer coefficient, g/sm 2 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

LWC Cloud liquid-water content, g/m 3 

MVD Water droplet median volume diameter, pm 


n Freezing fraction, dimensionless 

n 0 Stagnation freezing fraction, dimensionless 

p Pressure, Nt/m 2 

p w Vapor pressure of water in atmosphere, Nt/m 2 

p ww Vapor pressure of water at the icing surface, 

Nt/m 2 

r Recovery factor, dimensionless 

Re Reynolds number of model, dimensionless 

Res Reynolds number of water drop, dimensionless 

SLD Super-cooled large droplet 

s Distance along airfoil surface measured from 

stagnation line, cm 
tf Freezing temperature, °C 

t s Surface temperature, °C 

t Air temperature, °C 

T Absolute air temperature, K 

V Air velocity, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

We c Weber number based on model size and air 
properties, dimensionless 

We i Weber number based on model size and water 

properties, dimensionless 

Po Collection efficiency at stagnation line, dimen- 

sionless 

<j) Droplet energy transfer parameter, °C 

A Droplet range, m 

/^stakes Droplet range if Stokes Law applies, m 

Af Latent heat of freezing, cal/g 

A v Latent heat of condensation, cal/g 

p Air viscosity, g/m s 

6 Air energy transfer parameter, °C 

p Air density, g/m 3 

Pi Ice density, g/m 3 

p w Liquid water density, g/m 3 

a Surface tension of water over air, dyne/cm 

r Accretion time, min 

Subscripts 

st static 

tot total 


NAS A/CR— 2005-213853 


1 



Introduction 

Recent SLD-to-Appendix C scaling studies in the 
IRT 1 " have provided some important information as to 
how effectively existing scaling methods can be applied 
to scale SLD drop sizes to Appendix-C conditions. The 
test results showed that, within the IRT current capabil- 
ity in the SLD regime, good scaling can be achieved by 
matching scale and reference values of the parameters 
no and We L and the product /LA, provided model size 
ratios were limited to 3.4:1 and freezing fractions cov- 
ered the range from 0.3 to 0.5. Some of the results are 
reproduced in Figure 1. 

It was also found that main ice shapes appear to be in- 
dependent of MVD. Figure 2 shows ice shapes obtained 
by Chen' from a 61-cm-chord GLC 305 airfoil model 
test in the IRT in 1998 for MVD of 55 and 20 pm at 
freezing fractions of 0.3 and 0.5. At each freezing frac- 
tion, the values of A c , n 0 , (/>, 0 and Re were maintained 
constant as drop size decreased. Although //« changed 
from about 92 to 74%, there was no measurable effect 
of reducing MVD on the main ice shape. Similar results 
were observed later by the authors for NACA 0012 
airfoil model tests in the IRT. Figure 3(a) shows ice 
shapes obtained with 200 and 40 pm drop sizes at 51 
m/s for the 53-cm model, and 3(b) compares a MVD of 



-0.04 0.00 0.04 0.08 0.12 0.16 -0.04 0.00 0.04 0.08 0.12 0.16 

x/c x/c 

C = 91.4 cm. 1-25-02 Run 1 c = 91.4 cm. 2-8-02 Run 7 

c = 53.3 cm. 1 1-13-00 Run 1 c = 26.7 cm. 1 1-16-00 Run 6 

(a) Model Chord Scaled from 91 to 53 cm; n 0 , 0.3. (b) Model Chord Scaled from 91 to 27 cm; n 0 , 0.5. 
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0.49 

12.2 

14.2 
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Figure 1. Examples of Scaling from SLD Conditions to Appendix C. NACA 0012 Airfoils Tested in NASA Glenn 
Icing Research Tunnel at 0° AOA. 2 


200 pm ice shape with that of 20pm at 77 m/s for the 
91 -cm model. For each of the comparisons, the main 
ice shape was little changed by the reduction in MVD 
from SLD conditions to an Appendix-C value. These 
results suggest that it is possible within the range of test 
conditions given here to simulate SLD conditions with 
Appendix-C drop sizes, if the main ice shape is of pri- 
mary interest. 

However, two distinct features of those SLD ice shapes 
should be noted. First, when drop sizes were larger 
than 100 pm, the feather structures aft of the main ice 
shape were significantly larger than those of the Ap- 
pendix-C shapes. These large feathers were particularly 
prominent for a MVD of 200 pm. Second, the icing 
limit has changed with MVD (i.e. with J3 n ) because the 
SLD collection efficiency was so much larger than that 
for Appendix-C conditions. The SLD ice shapes fea- 
tured small feathers well aft of the Appendix-C icing 
limits. 

It is not easy to determine the exact conditions for large 
feather formation and the physics behind them in an 
icing tunnel environment. But this information is im- 
portant if one wants to know whether they can be simu- 
lated in small-droplet accretions. At rime conditions 
(i.e. n 0 = 1) the ice growth rate is mainly controlled by 
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Figure 2. Ice Shapes forMLD’s of 55 and 20 pm. 61-cm-Chord GLC 305 Airfoil Tested in NASA Glenn IRT at 0° 
AOA. 3 


the droplet impingement rate. This rapid freezing fea- 
ture at rime probably precludes splashing, and rime 
tests may therefore give some clues about the cause of 
the very large glaze SLD feathers. If we find that the 
large feather formations still form even in rime condi- 
tions, we would conclude either that splashing does 
occur in rime or that some other mechanism produces 
the large feathers. On the other hand, if the large feath- 
ers are not present in rime that would be strong evi- 
dence both that splashing is absent from rime and that 
splashing is probably the cause of the large feather fea- 
tures seen in glaze SLD accretions. 

Consequently, a study was commenced to record ice 
shapes produced at rime conditions with drop MVD’s 
ranging from a low of 15 pm to values well into the 
SLD regime. In the process of performing these tests, 
unexpected differences in the main ice shape became 
evident. In this paper, rime ice shapes obtained only for 
the MVD range of 15 to 100 pm will be presented to 
illustrate the differences in main ice shape observed. 
The results of tests with A/LD’s greater than 100 pm to 
explore large SLD feathers is beyond the scope of this 
paper, and plans are to present that information at an- 
other time. For the tests to be reported here, two 


NACA 0012 airfoil models with chords of 91.4 and 
53.3 cm were used at an airspeed of 77 m/s. 


Similarity Parameters 

The similarity parameters used in this study followed 
the Ruff 4 method. To be brief, only the final equations 
for the similarity parameters will be presented here. 
Interested readers are referred to a recent comprehen- 
sive review on icing scaling methods by Anderson and 
Tsao 5 for more details. 


The modified inertia parameter, K 0 , was defined by 
Langmuir and Blodgett: 6 




(1) 


In equation (1), K is the inertia parameter, 

K _ P« MVP 2 V 
18 d /j 


(2) 


For cylinders, the cylinder radius should be used in 
place of d in eq. (2). For airfoils, d is twice the leading- 
edge radius of curvature. For subsequent equations, d 
represents either the diameter for cylinders or twice the 
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(a) MVD, 200 and 40|im; c, 53.3 cm; V, 51 m/s. (b) MVD, 200 and 20|im; c, 91.4 cm; V, 77 m/s. 
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Figure 3. Comparisons of Ice Shapes from Appendix-C and 200-|im Sprays. NACA 0012 Airfoils Tested in 
NASA Glenn Icing Research Tunnel at 0° AOA. 2 


leading-edge radius for airfoils. For the NACA 0012 
airfoil model, a leading-edge radius of 0.0158c was 
used (see Abbott and von Doenhoff 7 ), where c is the 
model chord. MA Slo kes is the droplet range parameter, 
defined as the ratio of actual droplet range to that if 
Stokes drag law for solid-spheres applied. It is a func- 
tion only of the droplet Reynolds number, Reg. 

VMVDp ... 

Re S = ( 3 ) 

/' 


Langmuir and Blodgett’s tabulation of the range pa- 
rameter was fit to the following expression for this 
study: 


1 

^ Stokes 


'0.8388 + 0.001483 1 

v +0. 1847^ 


( 4 ) 


The accumulation parameter is: 


A, = 


LWCVt 

d Pj 


(6) 


If all the water impinging on the leading edge freezes at 
that location and the leading-edge collection efficiency 
is 100%, A c is a measure of the normalized thickness of 
ice that will accrete. 


The freezing fraction is defined as the ratio of the mass 
of water that freezes at a given location on the surface 
to the total mass of water that impinges the surface at 
that location. From Messinger’s 8 steady-state surface 
energy balance formulation, the freezing fraction is 


n = 


A r 


<!> + 



( 7 ) 


Of more practical interest than K 0 is the collection effi- 
ciency at the stagnation point, J3 0 , which was shown by 
Langmuir and Blodgett to be a function only of K 0 , 

X .S4 


1.40 K 0 


Po = - 


1 


1+1.401 K 0 - 1 - 


.84 


( 5 ) 


The individual terms in this expression are (/), the water 
energy transfer parameter, 

V 2 

V = t f -t st — ( 8 ) 

^ C p,ws 

9, the air energy transfer parameter, 
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6 = 


t , -c 


_ 

2c 


pJ 


Pww _ Ptot Pw 
T s t T tot p st 
1 Ptot Pww 


k (9) 


V-622 T to< 


L St J 


and h. the relative heat factor, introduced by Tribus, et. 
al. 9 


b = 


LWC V P 0 c p 

K 


(10) 


Equation (9) given by Ruff includes compressibility 
effects. Simpler forms without compressibility have 
also been used by Charpin and Fasso 10 and others, but 
the differences in values are not significant. For the 
present study the freezing fraction at stagnation, n (h will 
be used. At stagnation the recovery factor, r, in equa- 
tion (9) is unity. 

For completeness, the additional similarity parameters 
used in scaling studies will be given here. They are the 
Reynolds number for the model. Re, the droplet Weber 
number, We, the model Weber number using air den- 
sity, We c , and the Weber number based on model size 
and water density, We L 2 


RcV’‘' P 

P 

(11) 

We J lMVDp » 

(J 

(12) 

We c = V2dp 
a 

(13) 

V 2 dp 

We, = 

(14) 


a 


Test Description 


The models used were NACA 0012 airfoil sections with 
chords of 91.4, 53.3 and 26.7 cm. The 91.4-cm-chord 
airfoil is pictured in figure 4 (a). It was a full-span, 
fiberglass model. The 53.3-cm-chord model was of 61- 
cm span and made of aluminum. It was mounted verti- 
cally between splitter plates at the center of the IRT test 
section as shown in figure 4(b). Florizontal lines at the 
leading edge indicated tunnel center, ±2.5 cm and ±5 
cm from the center as visual guides for locating ice 
tracings. All tests were run at 0° AOA although the 
mounting system allowed rotation of the model for any 
given angle of attack setting. Also because of the quick 
start capability of the current IRT spray system, the 
models were not shielded during the initiation of the 
spray. 

In preparing for a test, the temperature and airspeed in 
the test section and the air and water pressures on the 
spray manifolds were set. When these conditions had 
stabilized, the spray nozzle valves were opened to initi- 
ate the spray. The spray was timed for the required 
duration, and then turned off. The fan was brought to a 
full stop and the researchers entered the test section to 
record the ice shape through hand tracings and photo- 
graphs. A thin heated stainless plate with a cutout in 
the shape of the airfoil leading-edge region was inserted 
into the ice to melt a thin slice down to the model sur- 
face. A cardboard template was placed into this gap 
and an outline of the ice shape traced. Tracings were 
taken at the vertical center of the tunnel (91 cm from 
the floor) and at 2.5 cm above the center. The tracings 
were digitized and the x-y coordinates for each ice 
shape recorded. Subsequently the coordinates were 
normalized by the model chord for comparison of ice 
shapes. 

Results from test entries in 2002 and 2003 will be pre- 
sented. Since the shape differences between the two 
tracing locations were never significant, only centerline 
shapes will be reported here. 


The icing tests were performed in the NASA Glenn 
Icing Research Tunnel (IRT). The IRT is a closed- 
loop, refrigerated, sea-level tunnel with a 1.8 m by 2.7 
m rectangular test section. It uses 10 spray bars, a con- 
figuration in operation since 1998, to generate a cloud 
of super-cooled water droplets. The Appendix-C cloud 
calibration used for these tests was performed in the 
summer of 2000. 11 

The SLD calibration used in this study was made in the 
summer of 2002 applying the same methods as the Ap- 
pendix C. At the time of the tests reported here, only a 
few specific MVD-LWC conditions for speeds of 51, 77 
and 103 m/s (i.e. 100, 150 and 200 kt) had been cali- 
brated. Therefore SLD tests were limited to these par- 
ticular conditions. 


Uncertainty Analysis 

Estimates of the uncertainty in the reported average 
conditions were made by considering inherent errors of 
instmments, temporal fluctuation and spatial variation 
of the instillment readings in the test section, and uncer- 
tainty in tunnel calibration of MVD and LWC. Total air 
temperature was believed to be good to ±0.5°C, and the 
uncertainty in air velocity was estimated to be ±lm/s. 
For Appendix-C conditions the net uncertainty in MVD 
was estimated at ±12%. For SLD conditions it may 
have been as much as ±20%. These uncertainties are 
not referenced to an absolute value of MVD, which is 
unknown. Repeatability and scatter in the LWC calibra- 
tion data suggests the uncertainty is about ±12% for 
both Appendix-C and SLD conditions. 
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(a) 91.4-cm-Chord NACA 0012 Model Installed in IRT 
Test Section. 



(b) 53.3-cm-Chord NACA 0012 Model Installed in IRT 


Test Section. 

Figure 4. Model Description. 

These uncertainties in the test parameters were used to 
estimate the following uncertainties in the similarity 


parameters for the Appendix-C tests: 9% in f} Q , 12% in 
A c , 13% in no, 3% in Re, 13% in We and 5% in We L . 
For the SLD tests the uncertainties were: 3% in /3 0 , 
12% in A c , 1 1% in n 0 , 3% in Re, 21% in We and 5% in 
We L . 

Results 

IRT Repeatability 

Rime icing tests were made over a number of MVD’ s, 
chord sizes and airspeeds while maintaining constant 
PnAc- They provide a good check on tunnel L WC cali- 
brations. Figure 5 compares ice shapes obtained at dif- 
ferent entry dates for MVD’s of 15 pm (5(a)) and 100 
pm (5(b)). Average recorded test conditions and corre- 
sponding similarity parameters for each of the tests are 
given in a table accompanying the figure. Tempera- 
tures were calculated to provide just rime icing condi- 
tions (i.e. the warmest temperature permitted for rime). 
Because the recorded conditions may be slightly differ- 
ent from the planned values for each test, the similarity 
parameters intended to be maintained may not always 
have matched exactly. The differences observed in the 
ice shapes were well within the uncertainty in LWC\ 
therefore, these results provide confidence that the IRT 
test conditions are repeatable for both SLD and Appen- 
dix-C conditions. 

The Change of Rime Ice Shape 

Rime ice shapes do not all have the same profile. The 
shapes shown in figure 5(a) for 15 pm sprays are no- 
ticeably different from those shown in figure 5(b) for 
100 pm. Both were obtained with the 91-cm-chord 
NACA 0012 model at a velocity of 77 m/s. Spray con- 
ditions for all the test results shown produce rime ice 
(no = 1) with the same leading-edge ice thickness (i.e., 
PoA c the same). Because of the different drop sizes, the 
values of for the two sets of test are significantly 
different: about 51% for the 15-pm case and 93% for 
100 pm. For a constant PqA c , this change in p 0 requires 
that the 15-pm tests be run at an A c nearly double that 
for the 100-pm. 

The 15-pm MVD cloud produces a spearhead-like rime 
ice shape shown in Figure 5(a). It features a small ice 
ridge of relatively compact ice at the center with adja- 
cent regions of long rime feathers growing into the flow 
direction to form relatively flat shoulders. Near the aft 
extent of those long feathers, a sudden transition to a 
region with much smaller and sparsely distributed 
feathers is observed extending back toward the icing 
limit. The shape and features were found to be repeat- 
able for different values of LWC, as shown in figure 
5(a). For the conditions of run 8 on 5-3-03, a tempera- 
ture of-23°C would have been cold enough to produce 
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MVD = 15pm. 05-03-03 Run 8 
MVD = 15tiin. 08-09-03 Run 3 


MVD = l()())iin. 10-10-02 Run 6 
A/I 'D = 100 pm, 03-02-02 Run 7 


(a) MFD, 15 pm; J3 0 , 51%. (b) MFD, 100pm; J3 0 , 93%. 


Date/Run 

c, 

cm 

°C 

V, 

m/s 

MFD, 

pm 

LWC, 

g/m 3 

T, 

min 

Po, 

% 

PtPc 

n n 

b 

$ 

°C 

0, 

°c 

Re, 

10 4 

We, 

10 3 

We c , 

10 3 

We L , 

10 6 

(a) 5-3-03/8 

91.4 

-26 

77 

15 

1.00 

20.3 

51.1 

1.82 

1.0 

0.50 

25.5 

32.4 

18.9 

1.37 

3.55 

2.65 

8-9-03/3 

91.4 

-21 

77 

15 

0.50 

40.8 

51.0 

1.82 

1.0 

0.25 

20.6 

26.9 

18.1 

1.36 

3.42 

2.61 

(b)l 0-1 0-02/6 

91.4 

-26 

77 

100 

0.68 

16.4 

93.2 

1.82 

1.0 

0.61 

25.6 

32.5 

19.1 

9.16 

3.59 

2.63 

3-2-02/7 

91.4 

-26 

77 

100 

0.68 

16.4 

93.2 

1.81 

1.0 

0.61 

25.7 

32.6 

19.1 

9.16 

3.59 

2.63 


Figure 5. Repeated Rime Ice Shapes forMFD’s of 15 and 100 pm. 91-cm-Chord NACA 0012 Airfoil at 0° AOA. 


rime, and other tests performed at this temperature gave 
results consistent with those of figure 5(a). The fact 
that the decrease in temperature to -26°C made no 
change in the appearance of the ice assures us that the 
shapes of figure 5(a) represent fully rime ice. Similar 
spearhead ice shapes have been observed and reported 
by Brunet and Guffond 12 on cylinders. 

In contrast, figure 5(b) shows a convex type rime ice 
shape for the 100-pm droplet MFD. The ice shape con- 
forms roughly to the leading-edge contour of the model 
with some small rime feathers (pointing into the flow 
direction) on the ice surface near the icing limits. 

Understanding the physics behind the different shapes 
produced for rime conditions is a necessary step toward 
insuring that scaling methods include the most signifi- 
cant factors affecting ice shape. Scaling methods de- 
veloped without a complete understanding may work 
well for tested conditions but be deficient for other 
situations not evaluated during method development. 
One step to understanding is to identify ranges of con- 
ditions for which different shapes occur. Additional 
rime ice shapes obtained from the same model and ve- 
locity as for figure 5 with intermediate MFD of 20, 30 
and 70 pm were examined and shown in figure 6. Two 
observations are noted. First, the convex type rime ice 
shape is maintained for MFD of 30 pm or larger. Sec- 
ond, there is a shape transition for MFD less than 30 
pm from ( 1 ) a convex type rime ice shape to ( 2 ) a 
wedge-like rime ice shape for a 20-pm MFD cloud. 


then to (3) a spearhead-like rime ice shape for a smaller 
MFD of 15 pm. In the transition, the large feather 
structures seem to become more important and relevant 
in developing the final ice shape. Close-up pictures of 
a 15-pm (spearhead type) and a 30-pm (convex type) 
rime ice shapes are provided in figures 7 and 8 for visu- 
alization purpose, and they show some span-wise varia- 
tions in features, possibly due to cloud non-uniformity 
effects that can not be fully captured by the 2D tracing. 

Finally, another series of rime ice shapes were exam- 
ined. These were obtained at the same velocity as for 
figures 5-8, but with a smaller model (53.3-cm chord). 
The results are given in figure 9. Due to the smaller 
model size used, the /F values for each drop size are 
higher than the corresponding values for figure 6. Now 
the 15-pm MFD (figure 9 (a)) no longer shows a spear- 
head shape, but rather a wedge shape very much like 
the 20-pm test of figures 6 (a) and 9 (a). Drop sizes of 
30 and 70 pm (figure 9 (b)), produced the same convex 
shapes observed for the larger model in figures 5 (b) 
and 6 (b). 

It is reasonable at first to speculate that the rime ice 
shape transition from the convex shapes of figures 5 
(b), 6 (b) and 9 (b) to the wedge (figure 9 (a) and the 
20-pm shape of 6 (a)), then to the spearheads (figure 5 
(a) and the 15-pm shape of 6 (a)) is a function of the 
stagnation collection efficiency Po and the correspond- 
ing P distribution on the model surface. If the change 
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MVD = 15(a.in. 05-02-03 Run 5 
MVD = 20(j.ni. 05-02-03 Run 7 


MVD = 30(oin. 08-30-03 Run 7 
MVD = 70(oin. 03-02-03 Run 6 


(a) Shapes for 

15- and 20-pm MVD Spray. 



(b) Shapes for 30- and 70-pm MVD Spray. 


Date/Run 

c, 

cm 

*st, 

°c 

v, 

m/s 

MVD, 

pm 

UVC, 

g/m 3 

T, 

min 

Po, 

% 

PiAc 

no 

b 

4, 

°C 

e, 

°c 

Re, 

10 4 

We, 

10 3 

We a 

10 3 

We u 

10 6 

(a) 5-2-03/5 

91.4 

-22 

77 

15 

1.00 

20.3 

51.0 

1.81 

1.0 

0.50 

21.8 

28.2 

18.3 

1.36 

3.44 

2.62 

5-2-03/7 

91.4 

-26 

77 

20 

1.00 

16.9 

61.3 

1.81 

1.0 

0.60 

25.1 

31.9 

18.7 

1.83 

3.51 

2.63 

(b) 8-30-03/7 

91.4 

-21 

77 

30 

0.65 

21.7 

73.2 

1.80 

1.0 

0.46 

20.6 

26.8 

18.3 

2.72 

3.48 

2.63 

3-2-02/6 

91.4 

-23 

77 

70 

0.60 

19.3 

89.6 

1.81 

1.0 

0.52 

22.7 

29.2 

18.8 

6.41 

3.55 

2.63 


Figure 6. Rime Ice Shapes fovMVD’s of 15 to 70 |om. 91.4-cm NACA 0012 Airfoil at 0° AOA. 


of rime ice shapes is correlated with the initial stagna- 
tion collection efficiency, p 0 , the present results would 
suggest that the rime ice shape is of convex type when 
Po > 0.8, wedge type when 0.6 < Po < 0.8, and spear- 
head type when p 0 < 0.6. 

To see how the P profile might differ for the 15- and 
100-|om shapes of figure 5, some calculations were 
made using LEWICE 2.2 13 for a clean 91.4-cm NACA 
0012 airfoil model. A simple monodisperse cloud was 
utilized. The numerical results, see figure 10, show that 
as the MVD decreases from 100 urn to 15 pm, the P 
curve is reduced in both magnitude and extent (im- 
pingement limits) on the model surface. The p 0 value 
has dropped almost by half (from 0.93 to 0.51) and the 
impingement limits have been reduced by a factor of 
four. 

The most striking difference in the two collection effi- 
ciency curves is the narrowness of the 1 5 -pm plot com- 
pared with that for 100 pm. This characteristic sug- 
gests that, in the absence of feathers and assuming the 
shape of the P curve does not change as ice accretes, ice 
would tend to form in only a very narrow region on the 
leading edge of the model, thus producing a spearhead. 
However, as ice accretes on the model, the resulting 
ice-shape change will influence both the shape of the P 
profile and the value of p 0 . Brunet and Guffond pre- 
sented photographs to illustrate how the shape of ice 
formations on cylinders changes with time as ice con- 


tinues to accrete. Brunet and Guffond’s ice shapes 
evolved from a convex shape similar to those of figures 
5 (b), 6 (b) and 9 (b), to a wedge-like shape like those 
in 9 (a) and the 20-pm shape in 6 (a), then to a spear- 
head ice shape closely resembling those of the 15-pm 
tests in figures 5 (a) and 6 (a). Although Brunet and 
Guffond did not provide conditions for these tests, the 
ice they show has the appearance of rime. 

The Brunet and Guffond results demonstrate the impor- 
tance of accretion time (or effectively the accumulation 
parameter A c ) on ice-shape characteristics. There is, of 
course, an interaction between the rime ice shape at any 
time and the collection efficiency. The changing nature 
of the ice shape creates a time-dependence in the 
collection efficiency, which in turn influences the 
evolving rime shape. As noted previously, the present 
tests were made with a constant value of the product 
PnA c , to insure that the leading-edge accretions would 
be the same size. Thus, as Po decreased, A c had to be 
increased to compensate. Consequently, it is difficult to 
separate the effects of collection efficiency and 
accumulation parameter on the developing ice shape. 
Additional testing is needed to record accretions with 
small droplets over a range of values of A c and also to 
look at larger droplets with greater values of A c than 
tested in this investigation. Such a study would help to 
determine how the shape transition process (i.e. con- 
vex-wedge-spearhead) correlates with A c and /A 
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(a) Side View. 


(b) Front View with Profile Tracing. 


(a) Side View. 


(b) Front View with Profile Tracing. 


Figure 7. Rime Ice on 91-cm-Chord NACA 0012 for MVD = 15 pm. Test of 5-2-03 Run 5. 
See Figure 6 (a) for Conditions. 


the MVD by itself 
can not completely 
represent the cloud 
dynamics and its 
influence on accre- 
tions, the actual 
drop size distribu- 
tion for different 
cloud MVD’s 
should be exam- 
ined. It was shown, 
for example in the 
IRT, that the 95% 
cumulative LWC 
droplet diameter 
(which is a good 
estimate of the 
maximum droplet 
diameter in the 
cloud) is about 23 
pm for a MVD of 
15 pm, but is 90 
pm for a MVD of 
20 pm and 260 pm 
for a MVD of 100 
pm. Thus, as MVD 
increases the drop- 
let size distribution 
becomes broader. 
It is possible that 
these distribution 
differences for 
various cloud 

MVD's have an 
effect on the rime 
ice shapes. If so, 
the effect on shape 
can be expected to 
apply to glaze con- 
ditions, as well as 
to rime. 

For glaze icing in 
SLD conditions, 
splashing may also 
have an important 
effect on ice shape. 
For conditions for 
which splashing is 
important, a new 

Figure 8. Rime Ice on 91-cm-Chord NACA 0012 for MVD = 30 pm. Test of 8-30-03 Run 7. way of representing 
See Figure 6 (b) for Conditions. the cloud to better 

characterize the 

Finally, there is another factor to be considered for the effects of the varying momenta of different drop sizes 
rime ice shape transition. That is the drop size distribu- ma y be needed, 
tion for different cloud MVD' s. Since it is known that 


NAS A/CR— 2005-213853 


9 





MVD = 15pm, 09-03-03 Run 4 

MVD = 20pm. 09-03-03 Run 5 

(a) Shapes from 15- and 20-pm MVD Spray. 


MVD = 30pm. 09-03-03 Run 8 

MVD = 70pm. 10-1 1-02 Run 7 

(b) Shapes from 30- and 70-pm MVD Spray. 


Date/Run 

c, 

cm 

ftt, 

°c 

V, 

m/s 

MVD, 

pm 

LWC, 

g/m 3 

T, 

min 

Po, 

% 

PcA-c 

n n 

b 

<f>, 

°C 

9, 

°C 

Re, 

10 4 

We, 

10 3 

We 0 

10 3 

We L , 

10 6 

(a) 9-3-03/4 

53.3 

-12 

77 

15 

0.50 

19.4 

62.7 

1.82 

1.0 

0.23 

11.6 

15.7 

10.0 

1.38 

1.96 

1.55 

9-3-03/5 

53.3 

-14 

77 

20 

0.50 

16.9 

71.6 

1.82 

1.0 

0.27 

13.0 

17.6 

10.1 

1.84 

1.96 

1.54 

(b) 9-3-03/8 

53.3 

-19 

77 

30 

0.65 

11.4 

81.4 

1.82 

1.0 

0.39 

18.1 

23.9 

10.4 

2.76 

2.01 

1.55 

10-11-02/7 

53.3 

-20 

77 

70 

0.60 

10.8 

93.1 

1.82 

1.0 

0.41 

18.9 

24.8 

10.7 

6.41 

2.04 

1.54 


Figure 9. Rime Ice Shapes for MVD' a of 15 to 70 pm. 53.3-cmNACA 0012 Airfoil at 0° AOA. 


Concluding Remarks 


Ice shape comparisons were made for NACA 0012 air- 
foil models with chords of 53.3 and 91.4 cm. Test 
conditions were set so that rime ice would be obtained 



Distance Front Stagnation Line, s/c 

Figure 10. Effect of MVD on Collection Efficiency, ft 
from LEWICE 13 for Clean 91-cm-Chord NACA 0012 
Airfoil at 0° AOA. t s „ -33°C; V, 77 m/s. 


ditions were set so that rime ice would be obtained for a 
nominal airspeed of 77 m/s, MVD ’ s ranging from 15 to 
100 pm and LWC’s of 0.5 to 1 g/m 3 . The product /3oA c 
was maintained the same for all tests (thus, the leading- 
edge ice thickness was the same) to permit direct com- 
parison of ice shapes obtained at different conditions. 
Important findings and suggestions are summarized as 
follows. 

• For the conditions of this study, a rime ice shape 
transition from convex to wedge to spearhead type 
was observed as the MVD decreased from 30 to 20 
to 15 pm. 

• For the conditions of this study, the convex type 
rime ice shape was maintained for MVD of 30 pm 
or larger. This suggests that it is possible in vari- 
ous rime icing conditions tested here to simulate 
SLD conditions with Appendix-C drop sizes no 
less than 30 pm. 

• For the conditions of this study, the rime ice shape 
was of convex type when ft 0 > 0.8, wedge type 
when 0.6 < Po < 0.8, and spearhead type when Po < 
0.6. 

• Evidence from reference 12 suggests that accretion 
time (or effectively the A c ) is important in deter- 
mining the ice shape. Further study of the effect of 
A c on the ice shape transition is needed. 
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• The actual drop size distribution for different cloud 
MVD’s may need to be considered to help explain 
the shape transition. 
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